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I. INTRODUCTION 



The electric dipole moments (EDMs) of electron and neutron are important observables 
for testing our ideas of CP- violation. In the standard model (SM), only one CP- violating 
phase exists in the Cabbibo-Kobayashi-Maskawa matrix. The predictions for the EDMs are 
extremely small, since the first nonzero contributions arise at two-loop level. They are several 
orders of magnitude smaller than the experimental limits Therefore, the EDMs are well 
suited for testing physics beyond the standard model 0]. 

In supersymmetric (SUSY) extensions of the standard model, additional CP-violating 
phases are possible. Moreover, the first nonzero contribution to the EDM already shows 
up at one-loop level. In particular, in the minimal supersymmetric standard model (MSSM) 
complex parameters can be introduced in the mixing matrices of squarks, sleptons, charginos 
and neutralinos, therefore yielding more possible sources of CP-violation. In weak-scale 
SUSY, the masses of the lightest supersymmetric particles are expected to be between 100 
GeV and 1 TeV. In this case, the EDMs can easily be much larger then the experimental 
limits and yield constraints on the CP- violating phases and on the other parameters of the 
MSSM. The conclusion would be that either the phases are small or the masses are large [0,0. 
Other arguments like the electroweak origin of the cosmological baryon asymmetry (BAU), 
however, would favor large CP- violating phases with relatively small masses [Q. Therefore, 
more careful analyses of the supersymmetric contributions to the EDM are necessary to 
clarify the situation. 

A suitable framework for numerical calculations in SUSY is the constrained MSSM, 
also called minimal supergravity-inspired model (mSUGRA) f^. In this model universality 
of the soft SUSY-breaking parameters at the grand unification (GUT) scale is assumed. 
Parameters of the model are the common gaugino mass M1/2, the common scalar mass Mq, 
the common trilinear scalar coupling parameter Aq, and tan/3 = V2/V1, with being the 
vacuum expectation values of the two Higgs fields. The number of independent complex 
phases can be reduced to two. The masses at the electroweak scale are determined by using 
renormalization group equations (RGEs). Such an approach to constrain the phases has 
recently been used in |^ and 

An important aspect in the calculation of the SUSY contributions to the EDMs of elec- 
tron (eEDM) and neutron (nEDM) is the fact that strong cancellations between the different 
contributions can occur. This has been particularly emphasized in , where the nEDM and 
eEDM in mSUGRA with two complex phases have been analyzed. Due to this cancellations, 
the bounds on the phases are less restrictive then those found in previous analyses. Addi- 
tional constraints on the phases originating from the cosmological bounds on the relic density 
of neutralinos have been studied in [|T^ . A different point of view has been presented in |TT[ 
where a model with seven independent phases at the electroweak scale has been assumed. 
Also in this analysis it was found that various cancellations between different contributions 
occur, and that in large regions of the parameter space the phases are not necessarily small. 

In our paper we analyze the eEDM and the nEDM simultaneously in mSUGRA with 
complex phases <y9^ and <^Ao 5 which are the phases of the higgsino mass parameter /i and the 
trilinear scalar coupling parameter Aq. We use RGEs to calculate particle masses, couplings, 
and phases at the electroweak scale from the input parameters at the GUT scale. We confirm 
the importance of cancellations. We find that quite general the cancellations occur between 
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the two most important contributions, which are the chargino and neutrahno contribution 
in the case of the eEDM and the chargino and gluino contribution in the case of the nEDM. 
Furthermore, the cancellations are only possible if the phases y?^ and ^a^^ are strongly 
correlated, in particular for small SUSY particle masses. In this case y?^ is strongly restricted. 
For the nEDM, there is also the problem of evaluating the hadronic matrix element. We use 
two different approaches, one based on the Quark-Parton Model 0], and a second one 



based on the Chiral Quark Model |]T3[. We find that the predictions for the nEDM are very 
different for the two models used. We show the regions in the M0-M1/2 plane which are 
excluded by the experimental bounds for both EDMs for specific values of the phases 9?^ 
and '^Aq- Finally, we also introduce an additional phase (y^s for the gluino mass parameter 
and study its influence. We find that also v^s is strongly restricted. 

In Sec. we give the expressions for the various contributions for eEDM and the quark 
EDMs, including the chromoelectric and purely gluonic dimension-six operator. We calculate 
the nEDM in terms of the quark EDMs in the two different models. In Sec. JTT|we determine 
the phases and MSSM parameters at the electroweak scale using the RGEs. In Sec. |V| we 
give the numerical analysis of the EDMs within mSUGRA and a discussion of the results. 
A summary is given in Sec. 0. Explicit forms of the mass matrices for sfermions, charginos, 
and neutralinos, as well as the expressions for the RGEs are given in the Appendices. 



II. CONTRIBUTIONS TO THE EDM OF ELECTRON AND NEUTRON 

The EDM of a spin-| particle is the coefficient of the effective operator 

He = -{^/2)dfh'a,JF^" . (2.1) 

We calculate the supersymmetric contributions to the EDMs of electron and quarks at one- 
loop level. In the case of the electron EDM we include chargino-sneutrino and neutralino- 
selectron loops. In the light quark case we include chargino-squark, neutralino-squark, and 
gluino-squark loops. For the chromoelectric dipole moments of quarks we include chargino- 
squark, neutralino-squark, and gluino-squark loops, whereas the gluonic dimension-six op- 
erator gets contributions from loops containing top quark, top squark, and gluino. 

The parts of the SUSY Lagrangian that are necessary to calculate the one-loop contri- 
butions mentioned above are 

^fxl U=9f i^kPn + bikPL) xl U , (2.2) 
^f'xtfm=9 f i^LkPR + ^LPl) Xt fm , (2.3) 
Cr,^^ = H9s/V2)qX'^ [e'^^^ni:,Pn-e-'^^^nZPL) ^ Qm , (2.4) 

where g and Qs are the electroweak and strong coupling constants, respectively, Pl,r = 
(1 =p 7^)/2, a = 1 . . .8 are the gluino color indices. A" are the Gell-Mann matrices, ip^ is 
the phase of the soft-breaking gluino mass. To simplify the notation the quark and squark 
color indices are suppressed. The scalar fields fi and Jr are linear combinations of the mass 
eigenstates fi^2- 
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/l \ 1-1 f f 



where TZ^ is the unitary diagonahzation matrix defined in Eq. ( |A8|) . Note that TZ^ depends 
on the phases ip^ and ipAo via the off diagonal entry of the squark mass matrix, see Eq. (|A5|) 
and Table |. The couplings are defined as (we use the notation of WW )'- 

= -5raiV,, , (2.6a) 
tj = -T^^iVn + YJli:,V,, , (2.6b) 
= -Ktu,, + Fe,d7^^^ , (2.6c) 



^mj ~ Ye^m,lUj2 i (2.7a) 

4, = , (2.7b) 

^■mj ~ YuTZ^iVj2 , (2.7d) 



^Lj — T^mlfij + '^tn2h^Rj ) (2.8a) 

4- = + ^^2/4- , (2.8b) 



hi- = r„(sin f5N^j - cos f3N^j) , (2.9a) 

hi^ = K(sin pn;^ - cos /?iv;.) = /i^;. , (2.9b) 

/^If = -V;,d(cos /?iV3,. + sin pN^j) , (2.9c) 

h% = -Fe,d(cos/3iV3*^. + sin/3iV;.) = /i^f , (2.9d) 



/4 = - [Q/ sin 2^H/A^rj + (1 - 2gj sin^ ew)N;^] /(v^cos 0^-) , (2.10a) 
fLj = [QfSm2ewN^j + (-2Q/ sin^ ^vy)A^2,] /(V2cos^h^) , (2.10b) 



v2mvK smp 

n,. = ^ ^"'^ ^ ■ (2.11b) 
V2mvi/ cosp 

Qf and are electric and Yukawa couplings of the fermion /, Oyy is the Weinberg angle, 
and tan/5 = V2/V1 is the ratio of the Higgs vacuum expectation values vi and V2- U and V 
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are the unitary matrices which diagonahze the chargino mass matrix, Eq. (|B1|). Naj is the 
unitary matrix which diagonahzes the complex symmetric neutrahno mass matrix, Eq. (|C^ ). 
For diagonahzing we use the singular value decomposition. 

A generic form for the one-loop EDM of spin-1/2 particles due to exchange of fermions 
and scalar particles has been worked out in Extensions of the EDMs to the full electric 



and weak dipole moment form factors for the top quark have been given in []16[. A non- 
vaninshing EDM demands a change in chirality of the external fermion and involves the 
imaginary parts of the couplings. In the following we give the complete analytic expressions 
for the individual one-loop contributions. We have compared our results with |^JTT| and 
found agreement. 



A. Chargino Contribution 

The chargino contribution to the EDM of the fermion / is given by 

where a = e^/(47r) and e = gsinOw- f is the isospin partner of / in the SU(2)-doublet. 
Neglecting the mass of the external fermions (in our case electron, up, down, and strange 



quark) the functions A and B have the simple form [15 



^W = ;T7T^('3-r + ^) , (2.13) 



2(1 -r)2 V 1 - r 

1 / 2r In r \ , , 

The first and second term in Eq. ( |2.12| ) are due to the Feynman diagrams Fig. la and Fig. lb, 
respectively. The expressions Im[r^] are given by: 

Im[r^^.] = lm[YMj2V,i]6rni , (2.15) 



Im [r:;,.] = Im [Y^V,2nUU,^nt - Y.Uj^nt)] (2.16) 
= (1/2)F„ ((1 - (-l)'"cos2^j)Im[f/,iK,2] +F,(-l)"^sin2^jIm[f/,-2V^,2e'^^^^ , 

id 

Im [F^J = Im [Y.U.^nUVjinZ - YuV.^nt:,)] (2.17) 
= (1/2)F, ((1 - (-I)'" cos2^ii) Im [Uj2Vji] + F„(-l)'" sm2e^ Im [f/j2V^,2e'^^]) . 
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B. Neutralino Contribution 



The neutralino contribution to the fermion EDM is given by 

4 2 , 



where 



k = l m=l Jm Jm 



-(l-(-ircos2^^;)Im[/^{ja 

~{l + {-ircos2ef)lm[h{ji,] . (2.19) 



C. Gluino Contribution 



The gluino contribution to the quark EDM is given by 

k=l Ik Ik 

2 2 
rv 777 ?77 

= -^sin((^3-¥^,-)sin20,-5^(-l)'=^g, B{^) , (2.20) 

k=l Ik Ik 

where = qI/^t^ and is the gluino mass. 

D. Quark Chromoelectric Dipole Moment and Gluonic Dimension Six Operator 

The quark chromoelectric dipole moment is defined as the coefficient dfl in the effective 
operator 

Cc = -{i/2)d'^qa^,^\X'l2)qG''>^' . (2.21) 

The chromoelectric dipole moment has also chargino, neutralino, and gluino contributions. 
They are given by 

4+ = - . E Im[ry^5(^) , (2.22) 

2 4 

m=l k=l y™ ™ 

and 
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where 



k = l <?fe 

= 3-flsin{^,-^^)sin2e,j2{-ir^C{^) , (2.24) 

k=l Ik Qk 



C(r) = 3A{r) - (l/3)5(r) . (2.25) 



The Wilson coefficient dc of the CP- violating gluonic dimension-six operator is defined 
through 

Cg = -{l/6)dGG^,aGl'G>;jabc . (2.26) 

The leading nontrivial contribution to dc in the MSSM is given by a two-loop diagram 
involving top, scalar top, and gluino [§,0: 



The definition of the two-loop function H can be found in [jT7|. 



E. EDM of electron and neutron 

Having defined the contributions from the individual Feynman diagrams, we can now 
write down the total EDM of the electron as the sum of neutralino and chargino contribu- 
tions: 

d' = dl++d\o . (2.28) 

In order to obtain the EDM of the neutron in terms of the quark EDMs, a specific 
description of the neutron as quark bound state is needed. Throughout this paper we use 
two different approaches. 

1.) The relativistic Quark-Parton Model: In this model, the contributions of the quarks 



to the nEDM are given in terms of quantities Ag |]12|, which are measured in polarized 
lepton-nucleon scattering: 

= r^^ (A^rf'^ + Adrf" + A,(i^) , (2.29) 

where the individual quark contributions are again given in terms of chargino, neutralino, 
and gluino contributions 

d" = 4+ + d\, + d] . (2.30) 

As already stated, the Ag are the measured contributions of the quark q to the spin of the 
proton; to use them for the neutron we have taken advantage of a simple isospin relation. For 
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definiteness we use the values given in Ref. |T^: A„ = 0.746, = —0.508, and = —0.226. 
The QCD correction factor t]^ takes into account that the quark EDM analysis is done at 
the electroweak scale and hence has to be evolved down to the hadronic scale with the help 
of RGEs. We use r]^ = 1.53 as given in Ref. |jl9 . 



2.) The Chiral Quark Model: This model is based on the effective chiral quark theory 
given in Ref. |]13|. The contribution of the quark EDMs to the nEDM is given by the 
nonrelativistic SU{6) coefficients 

r = (4/3)^'^ - . (2.31) 

The quark EDMs in this model are given by contributions of all quark and gluon operators 
(to leading order in as) with the proper dimensional rescaling. This yields 

= V^idl, + rf|o + dl) + V'^^idl. + d\, + 4) + ^^^c^G . (2.32) 

Tj^ , ^ and rf^ are the QCD correction factors due to RGEs, whereas Ag^ is the scale of 
chiral symmetry breaking in QCD; we use Vj^ = 1.53 [|19|, r^*^ ~ 7^<^ ~ 3.4 (as used in ||^), 
and AsB ^1-19 GeV 



III. DETERMINATION OF THE MSSM PARAMETERS AND PHASES 



The formulas for the EDMs, when evaluated in the MSSM with complex parameters in its 
most general form, contain too many free parameters. In order to study the constraints of the 
EDMs on the phases and mass parameters we have to reduce the number of free parameters 
by further theoretical assumptions. Therefore, we assume universality conditions for gaugino, 
sfermion, and Higgs mass parameters and the trilinear couplings 



Mo 

Ml/2 
^0 



Mi 



Ml 



= = Ml 
M2 = Ms , 
Ad, = 



Mr 



Mf 



(3.1) 
(3.2) 
(3.3) 



at the GUT scale Mqut 0, where i = 1,2,3 is the generation index. We determine the 
parameters at the electroweak scale with the help of the RGEs as given in . 

At the electroweak scale the following parameters can be complex: the trilinear couplings 
Af., the gaugino mass parameters M^., and the Higgs parameters fi and B. The product 
and the gaugino mass parameter M2 can be made real by redefinition of the fields. and 
B are determined by requiring the correct electroweak symmetry breaking: 



2/x5 



(ml, +ATi)-(m2,^ + Ar2)tan2/? 1 ^ 

K rriy 

tan^ /3 - 1 2 ^ 

(m|^ + m^^ + 2|/xp + ATi + AT2) sin 2/3 , 



(3.4) 
(3.5) 



where ATi 2 denote the leading one-loop corrections to the tadpole equations stemming from 
top, scalar top, bottom, and sbottom contributions p0|-p2|. The phase of fi, ip^, remains 
a free parameter. can be specified at any scale, because it does not evolve with the 
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corresponding RGE up to two loops . In order to determine the phases at the electroweak 
scale we assume M1/2 real, and Aq and n complex at the GUT scale. Note that at one-loop 
level only the phase difference between the phases of Aq and M1/2 is physically relevant. We 
summarize the complex phases entering the mass matrices in Table |. 

We use the following procedure for determining the soft SUSY-breaking parameters at 
the electroweak scale. We specify the gauge couplings, tan [3, and the Yukawa couplings of the 
third generation at the electroweak scale. We take Aq, Mq, M1/2 at Mqut with Eqs. ([3.11)- 
( p.3| ) as boundary conditions. The RGEs are given in the DR scheme. We evolve the RGEs 
for the gauge couplings at two-loop level from Q = mz to Q = Mqut which is determined 
by the condition Qi = g2- We evolve the RGEs for the Yukawa couplings at the one-loop 
level, because they enter the RGEs of the gauge couplings at two-loop level. We take into 
account threshold effects by including step functions for the coefficients of the beta functions 
(see e.g. [^]). For simplicity we assume that there is no mixing between the generations. 



We then evolve the RGEs for the soft SUSY-breaking parameters from Mqut to mz- The 
mass parameters Mj are decoupled from the RGEs if Mj{Q) = Q is satisfied. We calculate 
and B by requiring correct electroweak symmetry breaking Eqs. ( |3.4| ) and ( |3.5| ). The 



corrections are sensitive to the relative phases between the A parameters and /i. This phase 
dependence may change by a few GeV, which is in the range of the error expected by 



neglecting the other contributions to the one-loop corrected tadpoles pil , p2| . We iterate the 
complete procedure until the parameters vary less than 1%. 

For the discussion in the next Section it is convenient to have the following approxima- 
tions for the parameters at the electroweak scale at hand (the exact formulas for the one-loop 
results are given in Appendix 0). With acuT = 1/24 and Mqut = 2.38 x 10^^ GeV we get: 

M| ~ + 0.52M2/2 , (3.6a) 

M| ~ + 0.15M2/2 , (3.6b) 

M| ~ + 6.7Mi% , (3.6c) 

Mlc^Mlc::i Ml + 6.2M2/2 , (3.6d) 



Ml = M' ~ O.4IM1/2 , (3.7a) 
M2 = M ~ 0.82Mi/2 , (3.7b) 
M3 ~ 2.82Mi/2 , (3.7c) 



(l-l/)Ao-2Mi/2 , (3.8a) 

A^^ Ao- 2.8 Ml/2 , (3.8b) 

Ad-A^- 3.6 Ml/2 , (3.8c) 

Ae~Ao- 0.7 Ml/2 , (3.8d) 

where y varies between 0.85 and 1 for 40 > tan/3 > 1. Eqs. ( p.6a|) -( p.6d| ) are only valid 
for the first and second generation. Note that Eqs. ( p.8a] )-( p3.8d| ) have strong implications 
for the A parameters at the electroweak scale. If one takes, for example, M1/2 real and Aq 
imaginary, Aq = iA, at Mqut then one obtains the values given in Table 0. 
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IV. EDM ANALYSIS WITHIN MSUGRA 



In this Section we investigate the EDM of electron and neutron in the framework of 
mSUGRA with complex parameters. As outlined in Sec. |T|, this model is completely speci- 
fied by six parameters: Mq, M1/2, tan [3 and the phases ipAo and ip^. The experimental 
bounds obtained in [|| are \d^\ < dl^^ = 4.3 x 10"^'^ecm and \d^\ < d^^^ = 1.1 x 10~^^ecm. 

For the eEDM we have two supersymmetric contributions stemming from neutralino 
and chargino exchange, Figs, la and lb, respectively. The chargino contribution depends 
explicitly on the phase y?^, the dependence on (Paq comes only through the RGBs and is 
very weak. The neutralino contribution depends explicitly on and ipAo- the major 
part of the parameter space the chargino contribution dominates. The reasons are: (i) The 
loop function A{r), Eq. ( |2.13| ), entering in the chargino contribution, is larger than B[r), 
Eq. ( p.l4|) , which enters the neutralino contribution, (ii) The neutralino contribution is 
proportional to the selectron mixing angle sin26'e, which is usually rather small. 

In Fig. 2 we show d-o, the neutralino contribution of the eEDM, as a function of the 
CP-violating phases ip^ and (p^o with the other parameters fixed: Mq = 150 GeV, M1/2 = 
200 GeV, \Ao\ = 450 GeV, and tan/? = 3. As can be seen, the neutralino contribution alone 
already exceeds the experimental limit. The calculated eEDM is below the experimental limit 
only if cancellations between chargino and neutralino contributions occur. In this case the 
eEDM depends significantly on the phase ipAo if either Iv^/^l ■'^ \fAo \ot \Ae\ > tan p. In the 
first case the chargino contribution is small because it is proportional to sin y)^, therefore, the 
neutralino contribution can be of the same order of magnitude as the chargino contribution. 
In the second case the relevant phase in the neutralino contribution is determined by the 
off-diagonal element of the selectron mixing matrix Eq. (|X5|). In the mSUGRA model the 
absolute value of fi is fixed by the condition of radiative electroweak symmetry breaking 
Eq. ( |3.4| ). It turns out that \fi\ has always roughly the same order of magnitude as \Ae\ in 
the parameter region considered. Note that the neutralino contribution depends not only 
on the phase of {Ae — /i* tan/5), Eq. ( |A5| ), but also directly on via the neutralino mixing 
matrix, as can be seen in Eqs. ( p.l9| ), ( p.9c|) -( ^.llb|) . 

Due to the cancellation mechanism between chargino and neutralino contribution it is 
not straightforward to conclude which mSUGRA parameter values and phases are excluded 
by the experimental upper bound of the eEDM. To answer this question we show in Fig. 3 
the regions in the M0-M1/2 plane that are allowed by the experimental limit on the eEDM for 
different values of the phase 99^. In doing so we have taken ipAo = 7i"/2, which is the maximal 
phase difference between M1/2 and Aq at the GUT scale, tan/5 = 3, and |y4o| = 3Mq. For 
example, choosing ip^ = —0.1, the region in the Mq-Mi/2 plane to the left of the dashed- 
dotted line is excluded. As can be seen, the parameters Mq = 120 GeV and M1/2 = 160 GeV 
are allowed and give relatively light SUSY particle masses (for illustration: m^o = 58 GeV, 
m-± = 106 GeV, m^^ = 157 GeV, mg, = 160 GeV, mg, = 163 GeV). Taking = -0.54, 
only values of {M1/2, Mq) to the right of the solid line are allowed which, for example, 
means M1/2 > 0.9 TeV if Mq = 1.5 TeV or M1/2 > 1.5 TeV if Mq = 0.7 TeV. In this 
case rather heavy SUSY particles are predicted, (i. e. m^o > 398 GeV, m-± > 764 GeV, 
> 1170 GeV, rngj > 1073 GeV, mg^ > 1171 GeV). The bending in the dotted line 
for (ffj^ = —0.18 is caused by the cancellation mechanism between chargino and neutralino 
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contributions. The grey area is excluded, because the condition of radiative electroweak 
symmetry breaking is not fulfilled. 

Up to now we have only considered the eEDM. Now we consider the eEDM and the 
nEDM simultanously. Taking into account also the experimental upper limit on the nEDM 
will enlarge the excluded parameter region. The predicted value for the nEDM depends 
strongly on the neutron model which relates the nEDM to the EDM of its constituents. To 
demonstrate this fact we calculate the nEDM in the Quark-Parton Model and in the Chiral 
Quark Model as described in Sec. PE| . Also for the nEDM to fulfill the experimental bounds 
it is necessary that strong cancellations between the different contributions occur. 

Another way to show the systematics of these cancellations is to plot the allowed region 
in the ip^j,-ipAo plane. In the Figs. 4, 5, and 8 we consider rather small mSUGRA parame- 
ters: Mo = 150 GeV, M1/2 = 200 GeV, \Aq\ = 450 GeV, and tan/? = 3. In all iPf,-^Ao plots 
(Figs. 4, 7, and 8), the allowed values of the phases are within the small bands between 
the lines. In Figs. 4, 6, 7, and 8 we discuss the eEDM together with the nEDM. As can 
be seen from the dotted lines for the allowed region of the eEDM in Figs. 4a and 4b, is 
bounded, \ip^\ < 0.1, whereas ipAo is essentially unrestricted. However, the two phases have 
to be strongly correlated: for every (pAo, '^fj. can only vary in an interval Ayj^ < 0.01. Taking 
into account only the chargino contribution, one would obtain the restriction < 0.01. 

In Fig. 4a we show the experimentally allowed regions for the eEDM and the nEDM, 
calculated in the Quark-Parton Model. For the paramters chosen and the measured spin 



densities of the proton [|1^ A„ = 0.746, = —0.508, and A^ = —0.226, the allowed band 
in the ipf^-(pAo plane of the nEDM lies within the allowed band of the eEDM. In this case 
the nEDM is more restrictive. For the values of the spin densities taken, the nEDM and the 
eEDM have opposite signs. In Fig. 4b we plot the allowed band in the ip^-ipAo plane of the 
nEDM, calculated in the Chiral Quark Model, and compare it to the eEDM. They have the 
same sign. As one can see, in this case only a very small region of the parameter space is not 
excluded by experiment: \ip^\ < 0.01 and Iv^AqI ^ 0.15. (All phases have to be understood 
modulo TT.) 

In Figs. 5a and 5b we demonstrate the cancellation effects that play an essential role in the 
calculation of the nEDM. We choose the relation yj^ = — (tt/SO) ■ simpAo, which guarantees 
that the nEDM calculated in the Quark-Parton Model fulfills the experimental bound. We 
show the different contributions to the nEDM for the same parameters as in Fig. 4a. In 
Fig. 5a we show the corresponding chargino, neutralino, and gluino contributions. As can be 
seen, there is a strong cancellation between chargino and gluino contributions: each of the 
two contributions is approximately 18 times bigger than the whole nEDM. In Fig. 5b we show 
the up, down, and strange quark contributions to the nEDM. Again, cancellations between 
the individual quark contributions occur. It turns out, that the strange quark contribution 



is the most important one, as noted in |]12|. Therefore, it may turn out that an accurate 
measurement of the nEDM can also become a test of the spin structure of the neutron in 
the Quark-Parton Model. 

In the Chiral Quark Model the cancellations occur for up and down quark seperately. 
There are large cancellations between d'^+ and d^, between 0?^+ and rf^, and also between 
the resulting sums of this cancellations (see Eq. ( p.32| )). The purely gluonic dimension-six 
operator does not exceed the experimental limit by itself, however, it can further reduce the 
total nEDM. For the eEDM the cancellation between chargino and neutralino contribution 



11 



exhibits the same behaviour as shown in Fig. 5a, where the neutrahno contribution in the 
eEDM plays the same role as the gluino contribution in the nEDM. 

In Figs. 6a and 6b we show the regions in the M0-M1/2 plane which are excluded by 
simultanous consideration of the experimental limits on eEDM and nEDM. Fig. 6a is for the 
Quark-Parton Model and Fig. 6b is for the Chiral Quark Model. In both plots we choose 
the following values for the phases: (^9^0 = — vr/lO (dashed lines), y^yig = tt/S (dotted lines), 
ipAo = vr/2 (dashed-dotted lines), 93^ = — vr/lO (thin lines), and ip^ = — 7r/30 (thick lines). 
Only values of {M1/2, Mq) to the right of the corresponding lines are allowed by eEDM 
and nEDM simultanously. In Fig. 6a there are only four lines, because for = — tt/IO, 
ipAo = — vr/10 and ipfj_ = —tt/IO, = the parameter region Mq, M1/2 < 1-5 TeV is ex- 
cluded by experiment. The Quark-Parton Model is in general more restrictive than the 
Chiral Quark Model. However, in the Quark-Parton Model much smaller pairs of mass pa- 
rameters are allowed, for example Mq = 150 GeV and M1/2 = 200 GeV. As can be seen, the 
strongest cancellation effects are found for \(pAo | = vr/2 and sign pAo = — sign ip^. This is also 
observed in |jll|. If pAo and p^ have the same sign, the exclusion is more or less indepentent 
of Mo. 

Our numerical investigation of the nEDM includes the contributions of the one-loop 
gluino, chargino and neutralino exchange diagrams for the electric dipole operators. In the 
Chiral Quark Model we also include the chromoelectric dipole operators and the contri- 
bution of the purely gluonic dimension-six operator. In the following we want to discuss 
qualitatively which contributions are important to understand the behavior of the nEDM 
and its dependence on the mSUGRA parameters. 

The dominant contributions to the nEDM come from the chargino and gluino exchange 
diagrams of the quark EDMs. It is remarkable that the chargino contribution is almost 
independent of the phase pao- This is due to the fact that the second terms of Eqs. ( f^.l6[ ) 
and (|2.17| ) are suppressed by the Yukawa couplings Yu^d which are very small for light 
quarks. The gluino contribution, Eq. ( |2.20D , depends on both phases, and pAo: since it 
is proportional to the off-diagonal element of the squark mass matrix, rrig {Aq — /x*0(/9)), 
(see Eqs. (|A1|) and (|A5|) ) . The neutralino contributions to the quark EDMs are very small 
in contrast to the eEDM. 

In the Chiral Quark Model the down quark contribution is the most important one, 
because the EDM is proportional to {Add — du)- Moreover, for the chargino contribution 
we have Yd = (md/m„) tan /3 1^ > 6Yu if tan/3 > 3. The gluino contribution to the down 
quark is proportional to mdlTa[Ad — /i*tan/3], whereas the up quark EDM contains the 
factor rriu Im [A^ — /i* cot /?]. Taking into account that and \Aq\ have the same order of 
magnitude, we make the following observations: The down quark EDM depends mainly on 
/i. The up quark EDM is dominated by the term proportional to Au and is suppressed by a 
factor {iriu/md) cot P compared to the down quark term. 

The chromoelectric contributions (see Sec. |1I D| ) are suppressed by a factor gs/^An) com- 
pared to the electric dipole operator and, in general, they are less important. In the case 
where Mq > M1/2 the loop function C, Eq. ( |2.25| ), entering d?, Eq. ( |2.24| ), can compensate 
this suppression factor gg/^ATi). It also turns out that the contribution of the purely gluonic 
dimension-six operator is very small in the parameter region considered. 

In order to see how the restrictions on p^ and paq depend on the other parameters we 
also discuss a scenario with \Aq\ = Mq. In Fig. 7 we show regions in the p^-ipAo plane, 
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allowed by the experimental bounds on eEDM and nEDM in this case. We calculate the 
nEDM in the Quark-Parton Model with |^o| = = 150 GeV and the other parameters 
as in Fig. 4a. We find that the phase yj^o is less important than in the previous scenario 
{\Aq\ = 3Mq). The allowed values of 99^ are reduced roughly by a factor 1/3 compared to 
Fig. 4a, thereby suggesting a hnear dependence of the allowed values on l^ol keeping the 
other parameters fixed. Furthermore, the value of tan f3 effects the results in a similar way, 
because it enters in the off-diagonal element of the sfermion mixing matrix, Eqs. ( |A1D -( [A^ . 
This element is only important for the gluino contribution to the nEDM and the neutralino 
contribution to the eEDM. As can be seen in Fig. 7, the bands in the ip^-ipA^ plane, allowed 
by the eEDM and the nEDM in the Quark-Parton Model, overlap similarly as in Fig. 4a. 

In order to study the restrictions imposed by the universality conditions at the GUT 
scale, we modify the universality condition for the gaugino mass parameters, Eq. ( |3.2| ). We 
still assume M1/2 := Mi = M2 = IM3I, but introduce an additional phase (f^ for the mass 
parameter M3 at the GUT scale. We show in Fig. 8 the bands in the ip^j^-ifAo plane, allowed 
by the eEDM and the nEDM in the Quark-Parton Model, where we take for ^93 the values 0, 
7r/10, and 7r/5. We take the other parameters as in Fig. 4a. The eEDM depends on ^93 only 
via the RGEs, therefore, this dependence is very weak. Comparing the band of the eEDM 
(dotted line) with the bands of the nEDM for values of ^93 different from zero, one can see 
that (y93 is strongly restricted by experiment. A further possibility would be to introduce 
an additional phase ipi for the U{1) gaugino mass parameter Mi. This phase will enter the 
eEDM and the nEDM. It is expected that will change the restrictions on ^93 in a similar 
way as the phase ^9^0 changes the restrictions on ip^. 



V. SUMMARY 



We have studied the eEDM and the nEDM in the framework of mSUGRA with complex 
parameters. We have found that ip^ is strongly restricted by the experimental bounds. More- 
over, we have found that the phases ip^ and ipAo have to be strongly correlated, in particular 
for small values of the SUSY mass parameters, so that strong cancellations between the dif- 
ferent contributions occur. For the experimentally allowed values of the eEDM, the chargino 
contribution has to be cancelled by the neutralino contribution. The nEDM is dominated 
by the chargino and gluino contributions. The predictions for the nEDM depend very sensi- 
tively on the model which is used for the neutron. We have used the Quark-Parton Model 
and the Chiral Quark Model to calculate the nEDM. We have presented parameter regions 
in the M0-M1/2 plane which are excluded by simultanous consideration of the experimental 
bounds on the eEDM and the nEDM for different values of the phases <y9^ and ifAo- 
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APPENDIX A: SFERMION MASS MATRIX 



The sfermion mass matrices are given by 



l\/f^ — \ ILL JLR 

JLR JRR 



where 



Ml J + {Tf -Qf sin^ Ow) cos 2/5 m| + mj 



III 

Mj^^ = MJ^- + Qf sin^ Ow cos 2/5 m| + mj , 
Ml=mf\Aj-f^*Qm , 



Jlr 



with 



cot (3 for Tf 
tan P for Tf 



The eigenvalues are given by 

with ml < m\. We parametrize the mixing matrix TZ^ so that 

LA ^nf ( l"^^ ^( '^"^ V 



where yjj is given in Eq. (^^) and 



-Mj Mj - ml 
cose7 = < , sin^?=^^^^- >0 



Ilr' 



A^ = {Ml? + {ml-Miy 



APPENDIX B: CHARGING MASS MATRIX 



The chargino mass matrix 

j^x+ ^ [ ^ m,4/A/2sin/5 

cos (3 yU 

can be diagonahzed by the biunitary transformation 

where U and V are unitary matrices such that m~+ are positiv and m-+ 
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APPENDIX C: NEUTRALINO MASS MATRIX 



We define N^j as the unitary matrix which makes the complex symmetric neutrahno 
mass matrix diagonal with positiv diagonal elements: 



where m^o < m^o for ? < k. In the basis [24 



^a = {-tl,-iZ,H'',H'} , 
the complex symmetric neutralino mass matrix has the form 



/ m^y m 



\ 





uiaz mz 

mz /i sin 2/3 —/x cos 2/3 

\ -/icos2/? -/X sin 2/5 / 



where 



rriaz 



Mcos 

sin^VKCOs^w'(M-M') . 



(CI) 



(C2) 



(C3) 



(C4) 



APPENDIX D: SOLUTIONS OF THE ONE LOOP RGES 

The solutions of the one-loop RGEs (as given in pO[) for the soft SUSY-breaking pa- 
rameters are given by 

M.it) - ^ (Dl) 

MlW^MS^"-^lMt) (D2, 

Mlm^MU"-^[lMt)^lMt)) (D3) 

r2 



M|,«) = MS + [lut) + lut) + l^Mt)) (D6) 

, _ M,,, (hM - jmML) (07) 



AAt) 
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with 



Au{t) = ^ (^0 + Mt) - M„2H,{t)) (D8) 
A.[t) = - ^^^^ ,3(t) + 3Ut) + (D9) 

Ait) = A-^^^^(3Ut) + lM^^^ (DIO) 



t = HMgut/Q? (Dll) 

A = (D12) 

47r 

= 1^ (D15) 

16 3 13 

E{t) = {I + (3^t)—z {I + I32t) ^ (1 + (3it) W (D16) 

F{t) = [ E{s)ds (D17) 
Jo 

m) = ^ (jJsit) + 3Ut) + (D18) 
H2{t) = tE{t) - F{t) (D19) 
where hi = 11, 62 = 1, and 63 = —3. Some of the equations can be found in . 
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FIGURES 



(a) 



(b) 
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FIG. 1. The Feynman diagrams that contribute to the EDMs. 




FIG. 2. Ratio of the neutralino contribution d^o and the experimental hmit d^^^ of the electron 
EDM as a function of the phases ip^ and v^Aq- The mSUGRA parameters are Mq = 150 GeV, 
Ml/2 = 200 GeV, |^o| = 450 GeV, and tan/3 = 3. 
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FIG. 3. Boundaries of the areas in the Mq-Mi/2 plane excluded by the electron EDM, for the 
phases ipAo = vr/2, and ip^ = —0.54 (solid line), ip^ = —0.31 (dashed line), ip^ = —0.18 (dotted line), 
iff^ = —0.1 (dashed-dotted line). The areas to the left of the corresponding lines are excluded. The 
mSUGRA parameters are |^o| = 3Mo and tan/3 = 3. In the grey area the condition of radiative 
electroweak symmetry breaking is not fulfilled. 

Quark-Parton Model Chiral Quark Model 




FIG. 4. Bands in the (p^-(pAo plane allowed by the electron EDM (dotted line) and neutron EDM 
(solid line). The mSUGRA parameters are Mq = 150 GeV, M1/2 = 200 GeV, \Ao\ = 450 GeV, and 
tan P = 3. The neutron EDM is calculated in the Quark-Parton Model (a) and in the Chiral Quark 
Model (b). 
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FIG. 5. Cancellations of the various contributions to the neutron EDM in the Quark-Parton 
Model, taking the relation 93^ = — (vr/SO) • sin 93^0 . The mSUGRA parameters are Mq = 150 GeV, 
Ml /2 = 200 GeV, \ Aq\ = 450 GeV, and tan /3 = 3. (a) shows the chargino contribution (dashed 
line), neutralino contribution d^o (dotted line), gluino contribution (dashed-dotted line), and 
the whole neutron EDM (solid line), (b) shows the up quark contribution (dashed line), 
down quark contribution d'^ (dotted line), strange quark contribution d** (dashed-dotted line), and 
the whole neutron EDM (solid line). 
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Ml/2 [TeV] Ml/2 [TeV] 

FIG. 6. Boundaries of the areas in the M0-M1/2 plane allowed simultanously by the electron 
EDM and the neutron EDM. The neutron EDM is calculated in the Quark-Parton Model (a) and 
in the Chiral Quark Model (b). The phases are chosen as ipAo = — vr/10 (dashed lines), ipAo = it/5 
(dotted lines), ifAo = '^/'^ (dashed-dotted lines), and 99^ = — tt/IO (thin lines), 93^ = — vr/SO (thick 
lines). The areas to the left and below the corresponding lines are excluded. The mSUGRA pa- 
rameters are |Ao| = 3Mo and tan/? = 3. In the grey area the condition of radiative electroweak 
symmetry breaking is not fulfilled. In (a) the whole parameter region is excluded for 99^ = — vr/lO, 
VAo = -tt/IO and iff, = -vr/lO, 92^0 = 7r/5. 
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7r/50 




FIG. 7. Bands in the ip^-ipAo plane allowed by the electron EDM (dotted line) and 
neutron EDM (sohd line). The mSUGRA parameters are Mq = 150 GeV, M1/2 = 200 GeV, 
l^ol = ^0 = 150 GeV, and tan/3 = 3. The neutron EDM is calculated in the Quark-Parton Model. 
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FIG. 8. Bands in the ip^-ipA^ plane allowed by the electron EDM (dotted lines) and neutron 
EDM, calculated in the Quark-Parton Model, for ip^ = Q (solid lines), 993 = vr/lO (dashed lines), and 



ifs = Tr/5 (dashed-dotted lines). The mSUGRA parameters are Mq 
\Ao\ = 450 GeV, and tan/3 = 3. 
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TABLES 



TABLE I. Phases occuring in the mass matrices at the electroweak scale and at the GUT scale. 


mass matrix electroweak scale 




GUT scale 


Ml <^s = arg[^„-/x*cot/3] 
A/f? in arprf4j — /;*tan/?l 

M? ipe = BiglAe - n* tan /S] 






TABLE IL Values of the phases at the electroweak scale for M1/2 

at tho GUT scale when A = ;ril/j /2. 


real and Aq 


= lA imaginary 


X fAtly =0.85(1) 'PAu\y=0.85(l) 


"PAd 


'PAe 


0.1 -0.007 (0) -0.021 (-0.018) 
1 -0.075 (0) -0.203 (-0.177) 
10 -0.644 (0) -1.236 (-1.190) 


-0.0278 
-0.271 
-1.207 


-0.142 
-0.960 
-1.501 
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